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ABSTRACT 

Radio halos are synchrotron radio sources detected in some massive galaxy clusters. 
Their Mpc-size indicates that (re)acceleration processes are taking place in the host 
cluster. X-ray catalogues of galaxy clusters have been used in the past to search for ra¬ 
dio halos and to understand their connection with cluster-cluster mergers and with the 
thermal component of the intra-cluster medium. More recently, the Sunyaev-Zerdovich 
effect has been proven to be a better route to search for massive clusters in a wider 
redshift range. With the aim of discovering new radio halos and understanding their 
connection with cluster-cluster mergers, we have selected from the Planck Early source 
catalog the most massive clusters, and we have observed with the Giant Metrewave 
Radio Telescope at 323 MHz those objects for which deep observations were not avail¬ 
able. We have discovered new peculiar radio emission in three of the observed clusters 
finding: (i) a radio halo in the cluster RXCJ0949.8+1708; (ii) extended emission in 
Abell 1443 that we classify as a radio halo plus a radio relic, with a bright fila¬ 
ment embedded in the radio halo; (iii) low-power radio emission is found in CIZA 
J1938.3+5409 which is ten times below the radio - X-ray correlation, and represents 
the first direct detection of the radio emission in the “upper-limit” region of the radio 
- X-ray diagram. We discuss the properties of these new radio halos in the framework 
of theoretical models for the radio emission. 

Key words: Galaxy clusters; non-thermal emission; particle acceleration; radio emis¬ 
sion. Galaxy clusters: individual: Abell 1443, CIZAJ1938.3+5409, RXCJ0949.9+1708, 
RXGJ1354.6+7715 


1 INTRODUCTION 

Some galaxy clusters host diffuse radio emission which is not 
directly connected to any of the cluster radio galaxies. This 
emission fills the inner Mpc of the host cluster (radio halo) 
or is located at the cluster periphery (radio relics), and is 
characterised by a low surface brightness at z/ ~ 1.4 GHz (~ 
1 ^zJy/arcsec^ ), and steep radio spectrunfl with a >1 (see 
e.g. review bv iFeretti et al. li^. 

Radio halos are produced by ultra-relativistic electrons, with 
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^ We define the radio spectrum as S{iy) oc iy~^. 


Lorentz factors 7 l ~ 10^, spinning in large-scale fiG mag¬ 
netic fields. The radio power at 1.4 GHz (Pi.4GHz) corre- 
lates with the X-r ay luminosity (Lx) of the host cluster (e.g. 
Ibiang et al.ll2003 b Since Lx is a proxy of the cluster mass, 
this correlation indicates that more massive systems host 
more powerful radio halos. In additio n, a radio bi-modalit y 
is found in the Lx — Pi. 4 GHz plane (|Brunetti et al.|[2QQ9l b 
clusters with the same X-ray luminosity either host a radio 
halo, whose power follows the Lx — Pi. 4 GHz correlation (re¬ 
ferred as “radio-on state”), or do not host radio emission 
(“radio-off state”). 

From a statistical point of view, not only the power 
but also the occurrence of radio halos is known to increase 
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with Lx, and hence with the mass of th e host cluster (e.g. 
iGiovannini et al ] |l999l : ICuciti et al JIIH). 

T he physical mechanism powering this emi ssion is de¬ 
bated (|Keshet k, Loebll2010l : ICassano et al.ll2010h . However, 
a connection with the dynamical status of the cluster is es¬ 
tablished (e.g. iBuotel[ 2 ^ 001 ! : ICassano et aLllioiOl b 

Radio halos are likely formed during cluster mergers, 
that inject in the intra-cluster medium (ICM) a consid¬ 
erable amount of energy. Part of the energy is dissipated 
through turbulent motions and shocks. Turbulence could 
re-accelerate an existing population of seed r elativistic elec¬ 
trons and produce sy nchrotron emission fe.g. lBrunetti et ^ 
I2OOII: IPetrosianllioOll^ . 

Alternatively, it has been proposed that radio ha¬ 
los result from inelastic hadronic collisions between cos¬ 
mic ra y protons (CRp) and thermal protons (h adronic 
models, lDennisonlll980[ iBlasi &: Colafrances^ll999l b How¬ 
ever, hadronic models predict 7 -ray emis sion that has not 
been observed by the Fermi Satellite (e.g. lAckermann et al.] 
2014lb In some of t he hadronic model formulations (e.g. 
Keshet Loebll2010h they would require a different mag¬ 
netic held in clust ers with and without radio halos, which i s 
not observed fe.g. lBonafede et ^lOllalJCovoni et al.ll201Ql b 
In addition, hadronic models cannot e xplain radio halos wit h 
very steep spectral index {a >1.5, e.g. iBrunetti et al]l2008h . 


Mixed hadronic and r e-acceleration model s have 
been also propose d (e.g. IBrunetti &: LazarianI 1 20 111 . 
IZandanel et al.|[2014l b 

Although the pure hadronic models cannot explain the halos 
on the Lx — Pi. 4 GHzCorrelation, Mpc-scale diffuse emission 
originating from hadronic interactions is expected in the ra¬ 
dio off-state clusters too. The first evidence for it has been 
found by stacking ~ 100 images of massive clusters using 
the SUMMS survey ( Syd ney University Mol onglo Sky Sur- 
vev. iBock et al.l 1991) bviBrown et al.l (|201lh . According to 
IBrunetti Lazarian (|201lh . the electrons originating from 
hadronic collisions could produce radio emission that is a 
factor ~ ten below the Lx — Pi. 4 GHzCorrelation. The same 
electrons could be the seeding particles re-accelerated by 
turbulence during a merger. 

W e refer the reader to the review bv IBrunetti JonesI 
(|2014l i for a detailed description of the proposed mecha- 


Recently, radio halos have been found in clu sters with 
low X-ray luminosities (|Giovannini et ahl l 201 l| j. and we 
have discovered the hrst example of radio halo in a mas¬ 
sive coo l-core cluster that do es not show signs for a recent 
merger ([Bonafede et al.ll2014l b In addition, unexpected spec¬ 
tral features have been observed at low frequencies with 
early LOFAR observations, that have detected a steepen- 
ing of the radio spectru m at low frequencies in Abell2256 
([van Weeren et al]l 2012 l b Although the spectral measure is 
restricted to a small portion of the halo, it indicates that 
more complex models, either inhomogeneous turbulence or 
mixed hadronic and lept onic models, have to be considered 
([van Weeren et aLll2012l b 

Before the advent of Sunyaev-Zel’dovich (SZ) large sky 
surveys. X-ray catalogs were the primary source of infor¬ 
mation to “hunt” for new radio halos. However, due to the 
dependence of the X-ray luminosity on the square of the 
gas density. X-ray selected samples can be biased towards 



Figure 2. Zoom into the F-shaped source and relic region of 
A1443. The contours are the same as in Figin left panel. Colours 
refer to the WISE 3.4 /rm image. The colorbar is in magnitudes, 
not calibrated in terms of absolute surface brightness. 


low mass and cool-core systems. As noted bv iBas 3 (l2012t) . 
a better selection could be based on the SZ effect. 

The SZ signal integrated over the cluster angular 
extent measures the total thermal energy of the gas and 
as such correlates closely with the total cluster mass (see 
Planck collaboration, 2011, 2013 and references therein). 
Indeed, also the SZ signal c orrela t es str ongly with the 
rad io power, as disc o vered by Basul (l2012lj a nd co nhrmed 
by ISommer &: Basul (|2014| j: iGassano et al ] (l2013ll . After 
mergers, the thermalisation of the gas in the IGM happens 
on larger timescales than the boost in X-ray luminosity, 
while the SZ signal w ill grow moderately and gradually 
([Sommer fc Basul2014l b Hence, SZ surveys have the poten¬ 
tial of hnding radio halos in late mergers, and in clusters 
that are left out by X-ray selections. Recently, a radio 
bi-modality had been found in the Pi aghzSZ plane too , 
for clusters with M 500 ^ 5.5x lO^^M© ([Gassano et al.|[2013l b 


In this paper, we present a hrst sample of radio obser¬ 
vations conducted with the Giant Meterwave Radio Tele¬ 
scope with the aim of searching for new radio halos in 
clusters that have a strong SZ signal. We have selected 
from the Planck catalogues, pub l ished in 2011 and 2013 
([Planck Gollaboration et boill . l2013l j the most massive 
clusters (M500 > 6 x lO^^M©) for which no radio ob¬ 
servations were available to establish the presence of dif¬ 
fuse radio emission. In this paper we present the results 
from the hrst set of observations. The clusters presented in 
this work have been selected from the hrst Planck catalog 
([Planck Gollaboration et al.ll201lh . 

The remainder of the paper is the following: in Sec. [ 2 ] 
we describe the radio observations and the main steps of 
the data reduction. In Sec. [3] we analyse the results of the 
observations for each cluster. A discussion of the results is 
presented in Sec. 0] and we conclude in Sec. [5l Throughout 
the paper we use a AGDM cosmological model with Ho = 
71 km s“^ Mpc“\ Qm =0.27, Qa =0.73. 
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Figure 1. Abell 1443: Left panel: HR image at 323 MHz in colours and iso-brightness contours. The beam is 9.8”x8.1”. Contours start 
at 0.5 mjy/beam and are spaced by a factor 2. The first negative contour is dashed. The red dashed circle is centred on the X-ray centre 
of the cluster and has a radius of 1 Mpc. The sources subtracted in the LR image are marked by letters (A to G). Right panel: Colours: 
Exposure-corrected X-ray emission of the cluster from Chandra in the 0.5 - 7 keV energy range. The colorbar has units of photons cm“^ 
s“^ pixel“^. Radio iso-brightness contours from the radio LR image are overplotted. The beam is 27.4” x26.7”. Contours start at 1.1 
mJy/beam and are spaced by a factor 2. The first negative contour is dashed. 



Figure 4. CIZAJ1938.3+5409. Left panel: HR image at 323 MHz in colours and iso-brightness contours. The beam is 12.9” x8.8”. 
Contours start at 0.3 mJy/beam and are spaced by a factor 2. The first negative contour is dashed.The red dashed circle is centred on 
the X-ray centre of the cluster and has a radius of 1 Mpc. The sources subtracted in the LR image are marked by letters (A to H). Right 
panel: X-ray emission of the cluster from Chandra in colours and iso-brightness contours from the radio LR image. The colorbar has 
units of photons cm“^ s“^ pixel“^. The beam is 40.3” x34.9”. Contours start at 0.7 mJy/beam and are spaced by a factor 2. The first 
negative contour is dashed. 
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Figure 3. A1443: Grey colour scale: LR image of the cluster, 
blue contours are like in Fig. ^ right panel. Green contours are 
like in Fig. (T] left panel. Red contours are the residual of the HR 
image convolved with a beam as the one of the LR image. Red 
contours are drawn at ±3cr. The dashed-circle is centred on the 
cluster centre and has a radius of one Mpc. The restoring beams 
are shown in the bottom left corner. 
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Figure 5. GIZAJ1938.3+5409: Grey colour scale: LR image of 
the cluster, blue contours are like in Fig. [U right panel. Green 
contours are like in Fig. [U left panel. Red contours are the resid¬ 
ual of the HR image convolved with a beam as the one of the 
LR image. Red contours are drawn at ±3(7. The dashed-circle is 
centred on the cluster centre and has a radius of one Mpc. The 
restoring beams are shown in the bottom right corner. 


2 RADIO OBSERVATIONS AND DATA 

REDUCTION 

Observations were carried out at the Giant Meterwave 
Radio Telescope (GMRT) at 325 MHz, using a 33 MHz 
bandwidth subdivided into 256 channels and 8s integration 
time. Depending on the position of the target, the sources 
3G147, 3G48, or 3G286 were observed for 15 minutes at 
the beginning of the observing block, and u sed as absolute 
flux a nd bandpass calibrator, adopting the IScaife HealdI 
absolute flux scale. The absolute flux calibrators 
were also used to estimate the instrumental contribution 
to the antenna gains, which is needed for ionospheric 
calibration, and the instrumental phase information was 
used to correct the target field. In Table [2] further details 
about the observations are listed. 

In order to obtain a uniformly-reduced data sample, we 
have applied a semi-automated calibration procedure to the 
different datasets. This approach minimises the possibility 
of field-dependent effects introduced by a manual data re¬ 
duction. 

The main steps of the calibration procedure are out¬ 
lined below, and are based on AIPS (Astronom ical Image 
Processing Sy stem), SPAM (|lntema et al.ll2009l i and Obit 
(|Gottonl 120081 ) tools. Strong Radio Frequency Interferences 
(RFI) was removed from the target field data by statistical 
outlier flagging tools. Much of the remaining low-level RFI 
was modelled and subtracted from the data using Obit . Af¬ 
ter RFI removal, datasets have been averaged down to 24 
channels, in order to speed up the imaging process and, at 
the same time, avoid significant bandwidth smearing. 

For the phase calibration, we started from a model derived 
from the Northern Via Sky Survey (NVSS, iGondon et ahl 
1 9981) and WEsterbo rk Northern Sky Survey (WENSS, 
Rengelink et ahl Il997l ) when available. Using NVSS and 
WENSS, we have predcited the apparent flux of the sources 
at 323 MHz, using a simple spectral model fit and taking 
into account the GMRT primary beam. Then, a cut in flux 
is applied to keep only the stronget sources. Typically, 20-30 
sources are included in the first model. Then, we proceeded 
with self-calibration loops. We decided not to use a phase 
calibrator, as the GMRT field of view is wide, and a non- 
negligible flux is present in the field of the available phase 
calibrators. We note that using a phase calibrator to boot¬ 
strap the flux from the absolute flux calibrator would have 
altered the flux scale, leading to higher errors in the flux 
measurements. 

The first imaging steps have been performed using 
AIPS. In order to compensate for the non-complanarity of 
the array we used the wide-field imaging technique, decom¬ 
posing the GMRT field of view into ~ 80-100 smaller facets, 
depending on the field. We performed few rounds of cleaning 
and self-calibration, and inspected the residual visibilities for 
a more accurate removal of low-level REI. In order to cor¬ 
rect for ionospheric effects, leading to direction-dependent 
phase errors, we applied SPAM calibration and imaging to 
the target field. The presence of strong sources in the fleld- 
of-view enables one to derive directional-dependent gains for 
each of them (similarly to the peeling technique) and to use 
these gains to fit a phase-screen over the entire held of view. 
After ionospheric corrections, the last step of the calibration 
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procedure consists of subtracting all the sources outside the 
central 15^ from the target centre to facilitate the imaging 
steps. 

Processed data were imaged in CASA (the Common As¬ 
tronomy Software Applications package), and some more cy¬ 
cles of self-calibration have been performed. We have imaged 
the datasets at high resolution, using the Briggs weighting 
scheme (robust parameter equal to 0 in the CASA dehni- 
tion), and excluding the baselines shorter than 1.1 -0.9 kA, 
depending on the cluster’s redshift, to hlter out the possi¬ 
ble emission on Mpc-scale. High resolution images (HR im¬ 
ages) have been used to identify the radio radio galaxies in 
the cluster and to subtract their emission from the uv-data. 
Specihcally: we have selected the clean components corre¬ 
sponding to the radio galaxies detected in the HR images 
within 1 Mpc from the cluster centre, and we have sub¬ 
tracted the clean components from the visibilities. We have 
used these new datasets, where the sources’ emission had 
been subtracted, to search for diffuse radio emission in the 
clusters. To increase the brightness sensitivity to the dif¬ 
fuse emission, we have produced low resolution images (LR 
images) tapering down the baselines longer than 4 kA and 
using the Briggs weighting scheme (robust parameter equal 
to 0.3 in CASA definition). 

The final images of the full held-of-view centred on the 
target were corrected for the primary beam response. We 
estimate that the residual amplitude errors are of the order 
of 10% at 325 MHz in line with va lues reported for G MRT 
observation at this frequency (e.g. van Weeren et ah] 1 20091 : 
lintema et al.ll201ll : iBonafede et al.ll2012l b 


3 INDIVIDUAL CLUSTERS 
3.1 Abell 1443 

Little is known about Abell 1443 in the litera¬ 
ture. This cluster w as one of the first detec t ed b y 
the Planck mission |Planck Collaboration et al.l l201lh . 
IPlanck Collaboration et al.l ll201lh report a cluster mass, 
within rsocO of M 500 =7.7t°;® xIO^^Mq. The ROSAT 
Brightest Cluster Sample reports a n X-ray luminosity 
L[{)A- 2 AkeV] ~ 6.18 X 10^^ erg s“^ (jEbeling et al.l 120001 . 
corrected for the cosmological model used in this paper). 
Chandra archival observations are also available for this tar¬ 
get. We have retrieved the Chandra observations from the 
archive, and processed them in the standard way (see Sec. 
HI). The X-ray image is shown in Figure [1] The X-ray mor¬ 
phology appears elongated along East-West, and no dehned 
core is clearly visible. The gas distribution, as traced by 
the X-ray emission, is highly perturbed, indicating that the 
cluster is undergoing a merger event. 

As shown in Fig. [U we detect peculiar diffuse emission 
in the GMRT 323 MHz image. The cluster hosts a T-shaped 
source located at ~300 kpc from the cluster X-ray centre, 
which does not have the morphology of a radiogalaxy (Fig¬ 
ure o. Another patch of diffuse emission is located W of 
the T-shaped source (Figure [2|), located at ~700 kpc from 


^ rsoo is defined as the radius at which the cluster density is 500 
times the critical density. 
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Figure 7. RXCJ0949-hl708: Grey colour scale: LR image of the 
cluster, blue contours are like in Fig. [6] right panel. Green con¬ 
tours are like in Fig. [O] left panel. Red contours are the residual 
of the HR image convolved with a beam as the one of the LR im¬ 
age. Red contours are drawn at ±3f7. The dashed-circle is centred 
on the cluster centre and has a radius of one Mpc. The restoring 
beams are shown in the bottom left corner. 
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Figure 8. RXGJ0949.9+1708: colours and white contours refer 
to the emission at 1.4 GHz from the NVSS. White contours are 
drawn at (-1,1,2,4) mJy/beam. The beam is 45” x 45”. The first 
negative contour is dashed. Red contours refer to the LR GMRT 
image (as in Fig. [G] right panel). 


the cluster X-ray centre. After subtracting the sources A- 
G(Figure [T] left panel), we have imaged the data at low- 
resolution, as described in Sec. [2] 

Extended emission is detected, as shown in the low reso¬ 
lution image (Figure [U right panel). We note that no resid¬ 
ual emission from the sources A,B,C, and D is left in the 
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Figure 6. RXCJ0949.9+1708: Left panel: HR image at 323 MHz in colours and iso-brightness contours. The beam is 9.8” x7.9”. Contours 
start at 0.5 mJy/beam and are spaced by a factor 2. The first negative contour is dashed.The red dashed circle is centred on the X-ray 
centre of the cluster and has a radius of 1 M pc. The sources sub tracted in the LR image are marked by letters (A to H). Right panel: 
X-ray emission of the cluster from Chandra (|Ebeling et al.ll2010l) in colours and iso-brightness contours from the radio LR image. The 
colorbar has units of photons cm“^ s“^ pixel“^. The beam is 27.2” x26.8”. Contours start at 0.7 mJy/beam and are spaced by a factor 
2. The first negative contour is dashed. 
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Figure 9. RXCJ1354.6+7715. Left panel: HR image at 323 MHz in colours and iso-brightness contours. The beam is 19.4” x878”. 
Contours start at 0.5 mJy/beam and are spaced by a factor 2. The first negative contour is dashed.The red dashed circle is centred on 
the X-ray centre of the cluster and has a radius of 1 Mpc. The sources subtracted in the LR image are marked by letters (A to D). Right 
panel: X-ray emission of the cluster from Chandra 0.5 - 2 keV energy band in colours and iso-brightness contours from the radio LR 
image. The colorbar has units of photons cm“^ s“^ pixel“^. The beam is 33.7” x25.6”. Contours start at 0.8 mJy/beam and are spaced 
by a factor 2. The first negative contour is dashed. 
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Table 1. Radio observations. 


Planck name 

alt. Gluster name 

Obs. date 

RA and DEG 

z scale 

Flux cal 

O-HR 

O-RR 




h m s, d ^ 

kpc/" 

m 

Jy/beam mJy/beam 

PSZI G229.70+77.97 

Abell 1443 

27/01/2013 

12 01 27.7 +23 05 17.9 

0.27 4.I0I 

3C147 

0.13 

0.30 

PSZl G086.47+15.31 

CIZAJ1938.3+5409 

26/01/2013 

19 38 18.6 +54 09 33 

0.26 3.990 

3C48-3C286 

O.IO 

0.15 

PSZI G2I6.60+47.00 

RXCJ0949.8+1708 

25/01/2013 

09 49 51.7 +17 07 08 

0.38 5.175 

3C147 

0.15 

0.20 

PSZI GII8.46+39.31 

RXCJ1354.6+7715 

26/01/2013 

13 54 37.8 +77 15 35 

0.40 5.343 

3C286 

0.14 

0.30 

Col. 1: Planck name. Col. 

2: Other name from previous catalogues; Col. 3: Observation date 





Col. 4: Right Ascension and Declination of the target; Col. 5: target redshift; Col. 6: Source used to set the absolute flux scale; 



Col. 7: Angular to linear scale; Col. 8: rms noise in 

the high resolution image; Col 9: rms noise of the low-resolution image. 



Table 2. Radio sources 








Planck name 

Other name 

Detection 

Flux density Radio power at 323 MHz LAS x SAS 

LLS X SLS 




mJy 

WHz-i 

"x" 


kpc X kpc 

PSZI G229.70+77.97 

Abell 1443 

Whole emission 580+59 

1.4±0.1 Xl026 

330x230 


1350x940 



Halo(?) 

74.0 ±7.6 

1.7±0.2 xl025 

270x230 


1100x940 



Relic(?) 

57.0+5.8 

1.3±0.1 Xl025 

160x65 


660x270 



F-shaped source 450+45 

l.llto.l Xl 026 

100x60 


410x250 

PSZI G086.47+I5.3I 

GIZAJI938.3+5409 

Halo 

II.0+1.2 

2.4±0.3 Xl024 

180x100 


720x400 

PSZI G2I6.60+47.00 

RXGJ 0949.8+1708 

Halo 

21.0+2.2 

l.llto.l xlO^s 

200x80 


1040x400 


Col. 1: Planck name, Col. 2: Other name from previous catalogues; Col. 3: Extended source detected in the cluster; 

Col. 4: Flux density of the detected source, or upper limit to it. Col. 5: Radio power of the detected source assuming a spectral index ct = 1.3 for the k-correction; 
Col. 6: Largest and smallest angular scale of the source; Col. 7: Largest and smallest linear scale of the source in the adopted cosmological model. 



Figure 10. RXCJ1354.6+7715: Grey colour scale: LR image of 
the cluster, blue contours are like in Fig. [U right panel. Green 
contours are like in Fig. [U left panel. Red contours are the resid¬ 
ual of the HR image convolved with a beam as the one of the 
LR image. Red contours are drawn at ±3(7. The dashed-circle is 
centred on the cluster centre and has a radius of one Mpc. In the 
bottom left corner, the restoring beams are plotted. 


low resolution image, indicating that the emission from the 
sources embedded in the diffuse emission has been properly 
subtracted too. 

The Largest Linear Size (LLS) of the emission is ~ 1.3 
Mpc. It departs from the cluster centre toward W, encom¬ 
passing the F-shaped source and the patch of diffuse emis¬ 
sion, both detected in the high resolution image. Because of 


its peculiar morphology, a classification of this emission into 
the usual halo and relic categories is not straightforward. 
The Western patch of diffuse emission could be a relic, as 
it is at the edge of the X-ray emission and it has a LLS of 
~660 kpc. Given its morphology and the fact that no com¬ 
pone nt is visible in the WISE 3.4 image fFig.fTi lWright et m] 
l2010l i the possibility that we are seeing a radio galaxy is less 
convincing. 

The emission located centrally could be instead clas¬ 
sified as a radio halo, as it has a LLS of ~1.1 Mpc. How¬ 
ever, the emission itself is peculiar because of its morphology 
and it contains a b right F-shaped source. In another case 
(MA CSJ0717+3745, iBonafede et al.ll2009llvan Weeren et ahl 
l2009l i a brighter part, possibly a relic or a bright filament, 
has been found embedded into the radio halo emission. Spec¬ 
tral index and polarisation information are needed to dis¬ 
criminate the origin of the F-shaped source and understand 
the nature of the whole diffuse emission. Details about the 
radio properties of the emission detected in this cluster are 
listed in Tabled 

To verify that the diffuse emission is not affected by 
the incomplete subtraction of the sources detected at high- 
resolution (A-G in Fig.[T]) we have convolved the residual of 
the high-resolution image with a Gaussian having the major 
and minor axis as the restoring beam of the LR image (Fig. 
El). In this image, the HR contours are plotted over the LR 
ones, and over the contours of the residuals of the HR image 
convolved at low resolution. From this image we conclude 
that the sources A,B,G,D and F are properly subtracted. 
Some residual is detected corresponding to the position of 
source G, on top of the relic. This emission contributes with 
0.2 mJy to the relic emission. Therefore, it does not affect 
our estimate of the relic flux density. 
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3.2 CIZAJ1938.3+5409 

The cluster CIZAJ1938.3+5409 is located i n the so-called 
‘‘zone of avoidance”^ and was discovered bv lEbeling et al] 
(12003) by inspe cting X-ray data fr om the RASS Bright 
Source Catalog (|Ebeling et alJ [200G b The authors report 
an X-ray luminosity of T[o.i_ 2 . 4 fceV] ~ 10.89 x 10^^ 

erg s“^ (corrected for the cosmological model adopted 
i n this work). The cluster mass, as reported by Planck 
(|Planck Collaboration et ahl l2013h is M 500 = 7.51q 3 x 
lO^^M©. Therefore, it is quite a massive cluster, but un¬ 
fortunately nothing is known in the literature about its dy¬ 
namical state. 

Chandra observations are available for this target (Jones 
et ah, 2015, in preparation). The X-ray and radio emission 
of the cluster are shown in Fig. [U 

After the source subtraction, we discover extended ra¬ 
dio emission, located at the cluster centre and having a 
size of ~ 720 kpc (Fig. [H left panel). The flux density of 
the emission is 11 mJy, which translates into a power of 
2.4 ± 0.3 X 10^^ W Hz-^ at 323 MHz. The size and the 
position of the emission would classify it as a radio halo. 
We note that the corresponding power at 1.4 GHz would 
be Pi. 4 GHz ~ 4 X 10^^ W Hz“^, which is unusually low for 
such a massive clustefl Indeed, the halo is at least 8-10 
times under-luminous in radio for its SZ signal and X-ray 
power, respectively (see Fig. HU). The radio power at 1.4 
GHz would place the halo in the “upper-limit” region of the 
Pi. 4 GHz - M500 plane. 

To verify that the diffuse emission is not due to the 
the incomplete subtraction of the sources detected at high- 
resolution (A-H in Fig.[4|) we have convolved the residual of 
the high-resolution image with a Gaussian having the ma¬ 
jor and minor axis as the restoring beam of the LR image 
(Fig. [ 5 ]). In this image, the HR contours are plotted over the 
LR ones, and over the contours of the residuals of the HR 
image convolved at low resolution. From this image we con¬ 
clude that the sources A,B,G,D,F,G and H are properly sub¬ 
tracted, as no residual emission is detected neither in the LR 
image nor in the residual image. Some residual emission is 
detected corresponding to source E. This emission accounts 
for 0.2 mJy and it is only confined to the source position, 
hence it cannot explain the 11 mJy emission on larger scale, 
detected in the LR image. 


3.3 RXCJ 0949.8+1708 


The cluster RXGJ 0949.8+1708, also known as MAGS 
J0949.8+1708, is an X-ray luminou s cluster, discovere d 
in X-rays by the ROSAT satellite (|Ebeling et al.1 I 2 OOOI L 
The Planc k satellite detected the c l uster through the 
SZ effect (jPlanck Gollaboration et I 2 OIII L and found 
it to be massive, with M 500 — 8.2 ± 0.6 x lO^^M© 

(|Planck Gollaboration et al.ll20l3 L 

Ebeling et al ] (l20ld) have analysed the cluster in detail 
using the UH 2.2m optical telescope and X-ray Chandra ob¬ 
servations. They derive an X-ray luminosity Lx[o.i- 2 AkeV] — 
10.6 ± 0.6 X lO^^ergs”^ within P 500 , and a temperature 
T = 8.9 ±1.8 keV. 
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Figure 11. The morphological diagrams c — w (top), c — PsjPo 
(middle) and w — P — 3 /Pq, adapted from Cassano et al. (2010). 
Red filled symbols refer to radio halos (red filled circles are from 
Cassano et al. (2010), red empty circle is CL1821+643 from 
Bonafede et al (2014), red filled stars are the radio halos dis¬ 
covered in this paper. Blue filled dots refer to to mini halos, black 
circles are clusters with no radio emission (from Cassano et al. 
2010). The green star and arrows refer to the low power radio 
halo in CIZAJ1938.3+5409. Arrows are upper limits to the P3/P0 
value. 
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Figure 12. Correlation between the radio power of radio halos 
at 1.4 GHz and the cluster mass M500 as derived by SZ mea¬ 
surements (|Planck Collaboration et al.ll20l3 b Blue points are the 
radio halos used to derive the best-fit line (shaded line). Green 
points are halos with spectral index a > 1.5. Arrows are upper 
limits (points, arrows and bet-fit lines taken by see Cassano et al. 
2013). Red filled stars are the radio halos presented in this paper, 
for which a spectral index o; = 1.2 is assumed. The red arrow 
refers to the upper limit on RXCJ1354.6+7715 (halo size ~ 700 
kpc, SsigmaLR detection threshold, see text for details). 


In Fig. [6] the X-r ay emission of the c luster is shown. 
The MACS clusters in lEbeling et al.l (|2010l ) have been clas¬ 
sified with a morphological code from 1 o 4, 1, meaning pro¬ 
nounced cool-core and good optical-X ray alignment, and 
4 meaning multiple X-ray peaks and no cD galaxy. RXCJ 
0949.8+1708 is given a morphological code 2, meaning that 
although no pronounced cool-core is present, the X-ray emis¬ 
sion shows a good optical-X-ray alignment and concentric 
X-ray contours. 

_ This cluster has already been observed bv I Venturi et al.l 

(|2QQ8h at 610 MHz (the author refer to the cluster with the 
name Z2261), who found positive residuals after the subtrac¬ 
tion of the sources and indicated the cluster as a candidate 
radio halo. 

In Fig. [6] the radio image from our new GMRT obser¬ 
vations is shown. We detect diffuse radio emission located 
at the cluster centre and with a LLS o f ^1 Mpc, conf i rming 
that the positive residuals detected bv IVenturi et al.l (|2008l ) 
are part of more extended radio emission. We classify the 
emission as a radio halo. It is elongated in the SW-NE di¬ 
rection, and it does not follow the emission of the gas. 

To verify that the diffuse emission is not affected by 
the incomplete subtraction of the sources detected at high- 
resolution (A-H in Eig. [G]) we have convolved the residual 
of the high-resolution image with a Gaussian having the 
major and minor axis as the restoring beam of the LR image 
(Eig. [7|). Erom this image we conclude that the sources are 
properly subtracted. Some residual is detected on top of the 
brightest region of the halo. Likely, it is the brightest part 
of the radio halo, which is on a scale smaller than 1 Mpc 


and, as such, is not filtered-out in the HR image. The flux 
density corresponding to this emission is ~ 0.5 mJy, hence, 
even if it came from the individual radio sources it would 
not affect the estimate of the halo flux density (see Tab. [2]). 


3.4 RXCJ1354.6+7715 

RXGJ1354.6+7715, also known as MAGS J1354.6+7715, 
is an X-ray l umino us galaxy cluster discovered by 
iBohringer et al.l (|20Q0I L They report an X-ray luminosity 
in the energy band 0 .1-2.4 keV of Lx[o.i- 2 AkeV] — 9.4 x 
10 ^^ergs“^ (corrected for the cosmological model used in 
this work). 

The Planck satellite detected the cluster 
throu gh the SZ effect dPlanck Gollaboration et al.l 
l 201 lh . IPlanck Gollaboration et al.l ( 20131 ) report 
M 50 0 = 6.2 + 0.6 X lO^M© . 

iHoresh et al.l (|2010l ) have analysed the matter sub¬ 
structure in the cluster using strongly lensed arcs detected 
through HST observations. Their analysis suggests the exis¬ 
tence of two separate galaxy concentrations. They conclude 
that this cluster could be during some stage of a merger, 
with a considerable amount of substructure. 

Chandra observations are available in the archive. We 
have calibrated them in the standard way (see Sec. and 
in Eig.[9]the X-ray emission is sh o wn. In agreement with the 
optical analysis bv iHoresh et al.l (|2010l ). two gas concentra¬ 
tions are clearly detected. 

In Eig.[9]the radio emission from the cluster at 323 MHz 
is shown. A bright radio galaxy (G in Eig. (9] left panel) is lo¬ 
cated at the centre of the cluster. Its total flux is ~159 mJy. 
It has a nucleus and a lower surface brightness tail. Residual 
calibration errors are present around the radiogalaxy. The 
presence of the tail and the residual calibration errors make 
the subtraction of the clean components of the radiogalaxy 
G difficult. 

Some residual emission is present in the low-resolution 
radio image, which is likely due to the convolution at lower 
resolution of the residual calibration errors of the radio 
galaxy G and to residual emission from the tail. In order 
to confirm this, we have convolved the residuals, obtained 
after the cleaning of the high-resolution image, with a Gaus¬ 
sian having the major and minor axis as the restoring beam 
of the low-resolution image. The comparison is shown in Eig. 
M Hence, despite the massive and likely merging cluster, no 
radio emission in form of radio halo or relic is present at the 
sensitivity level reached by our observations. 

It must be noted that RXGJ1354.6+7715 is the least 
massive clusters among those analysed in this paper. Putting 
an upper limits to the possible radio emission below our 
detection threshold is not trivial, because of the presence 
residual calibration errors connected with the radio galaxy 
G, and more generally because of the non-universal surface 
brightness distribution of radio halos. However, considering 
a conservative detection threshold of 3a lr over a circle of 
350 kpc radius, we can exclude the presence of radio emission 
on a scaled 700 kpc with a flux density S ^ 5 mJy, which 
translates into an upper limit of Pi. 4 GHz < 5 X 1023. This 
value is a factor ~5 below the Pi,4ghz-SZ correlation. 
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4 DISCUSSION 

4.1 Substructure analysis 


All the clusters analysed in this work have Chandra archival 
observations. X-ray images are a powerful tool to understand 
the cluster dynamical state, complementary to optical ob¬ 
servations. 

We have reprocessed the Chandra observation of the 
clusters using the latest calibration tab les (CALDB 4.5.9) in 
Cl AO. Following ICassano et al.l (|2010l b the event files have 
been binned by a factor four in order to undersample the 
Chandra PSF. We have selected the 0.5-2 keV energy band, 
and we have verified that more than 2000 counts were de¬ 
tected inside the region of interest (a circle of 500 kpc radius) 
to allow for a morphological analysis. The images have been 
normalised by the exposure map, and the sources in the clus¬ 
ter field have been visually identified and excluded from the 
analysis. 

In order to quantify the dynamical status of the clusters, 
we compute three dynamical indicators, e xtensively used in 
the literature (e.g. iBohringer et al.ll2010h : the power ratio 
P3/P0, the centroid shift, w, and the concentration param- 
eter^_a _ 

ICassano et al.l (|2010h have used these methods, hnding 
that clusters with and without radio halos occupy different 
regions in the morph ological diagrams. Re cently, we have 
found a first outlier ([Bonafede et al.l l2014l i in the case of 
CL1821+643. 


The morphological indicators P3/P0, le, and c have been 
computed within an aperture radius Rap = 500 kpc, and 
considering as centre the centroid of the X-ray emission 
(P3/P0 computation) and the cluster X-ray peak (c and w 
computatio n), to have a fair com parison with the clusters 
analysed bv ICassano et al.) (|2010l i. We note that the clusters 
presented in this work ha ve mass and redshift comparable 
to the clusters analysed bv ICassano et al ] (|2010l i , hence the 
comparison makes sense. The errors and the photon bias 
have been estimated t hrough Monte-Carlo sim ulations, fol¬ 
lowing the approach of IBohringer et ah (l2010h. _ 

The centroid shift method ([Maughan et al.ll2008l i mea¬ 
sures the standard deviation of the projected separation be¬ 
tween the X-ray peak and the centroid in units of Rap, com¬ 
puted within circles of increasing radius (from 0.05x Pap to 
Rap, in steps of 0.05xPap). The centroid shift, re, has been 
derived computing the X-ray weighted centroid within cir¬ 
cles of increasing radius. 

The concentration parameter, c, measures the ratio be¬ 
tween the X-ray surface brightness within 100 kpc ov er the 


surface brightness within 500 kpc (jSantos et al.l 120081 1. It is 
higher for cool core and regular clusters, while it is lower for 
merging clusters. 

The power ratio is a multipole decomp osition of the 
two-d imensional projected mass distribution (jBuote &: Tsail 
Il995l b We have determined the power ratio P 3 /P 0 , which 
is the lowest power ratio mo ment providing a clea r sub¬ 
structure measure. Following IBohringer et al ] (l20ldl , mo¬ 
menta are computed excising the inner 0.05xPap. In the 
computation of P 3 /P 0 , the photon bias takes into account 
the spurious amount of substr ucture that could cont ribute 
to the power ratio. As noted bv IBohringer et ^ (|2010l b also 
very regular clusters would show some substructure because 
of the photon noise. Since we are comparing here with the 


analysis by ICassano et al.l (|2010l i , the bias will not be taken 
into account in the following analysis. However, the photon 
bias gives a larger contribution to regular clusters than to 
perturbed ones. Hence, once it is taken into account in the 
power ratio analysis, the separation of merging and regular 
clusters along the P3/P0 axis would likely be more evident. 

In Fig. [11] we show the position of the clusters in the 
P3/P0 — c, and w — c d iagram s, taken from ICassano et al.l 
(l2010l) . lBonafede et alJ (l2014h . 

4.1.1 Abell 1443 

The substructure analysis confirms the visual impression 
from the X-ray image. We obtain w = 0.0406 ± 0.004, 
c = 0.112 ± 0.006, and P3/P0 = (1-6 ± 0.7) x 10“®. The 
dynamical status of the cluster is highly perturbed, and the 
cluster falls in the merging quadrant of the P3 — PO, w, c di- 
agramfl. As expected for such massive and merging clusters, 
diffuse radio emission has been discovered. 


4.1.2 CIZAJ 1938 .3+5409 

We obtain an upper limit for the value of P3/Po < 3.8x10“®, 
and a rather high value of c = 0.292 ± 0.006 that are typ¬ 
ical for non-merging clusters. Instead, the shift parameter 
w = 0.040 ± 0.003 is typical of merging clusters. Although 
a major merger can be reasonably excluded, the analysis is 
not conclusive of the cluster dynamical status. 

In Fig. [11] the cluster position in the morphological 
planes is shown. We note the peculiar position of the cluster 
in the diagrams. While in the plane c — P3/P0 the clusters 
sits among the non-merger ones, it lies in an almost empty 
quadrant in the c — w and w — P 3 /PQ planes. 

We note that the cluster morphologica l indicators are 
very similar to the one of CL1821+643 ([Bonafede et al.l 
l2014l i. However, the radio emission in the two clus¬ 
ters are rather different. CL1821+643 hosts a radio 
halo with a power that follows the Pi. 4 GHz — Lx and 
Pi .4 GHz — SZ correlations, while the radio halo discov¬ 
ered in CIZAJ1938.3+5409 is at least 8 times below the 
Pi .4 GHz — SZ correlations. 

Recently. iKhatril ([2015l i has shown that Planck SZ mea¬ 
surements might be contaminated by CO Galactic emis¬ 
sion. This would result in an over estimate of the clus¬ 
ter mass from SZ measurements, that can be relevant 
for clusters close to the Galactic plane. In the case of 
CIZAJ1938.3+5409, a mass M500 ~ 4x lO^^M© would bring 
the radio halo on the Pi.4GHz — SZ correlation. As the error 
on Planck masses are not quantified yet, a more precise esti¬ 
mate of the possible error in this cluster cannot be derived. 
However, we note that the cluster doe s not fa ll in th e most 
CO contaminated regions derived bv iKhatril ([2015l b This 
suggests that the mass derived by Planck is a good estimate 
of the cluster mass. 

Since the c luster total Lumino sity in the band 0.1-2.4 
keV is known ([Ebeling et al.l[ 200 I . see Sec. 13.2D . we have 
checked the position of the radio halo in the Piaghz — Lx 


As expected, in this case the bias would not change the value 
of P 3 /P 0 dramatically. Taking the photon bias into account one 
would obtain P 3 /P 0 = (1.5 ± 0.7) x 
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plane (|Kale et al.ll2015h . Also in this diagram the radio halo 
is a factor ~ 10 below the Piaghz — Lx correlation, and 
falls in the upper-limit region. This provides an independent 
evidence - with respect to the Planck mass measurement - 
that the cluster hosts a radio halo of low power for its mass. 

In the next Section, we discuss the possible origin for 
the low power radio halo in CIZAJ1938.3+5409. 


4.1.3 RXCJ 0949 . 9 -hl 708 

For RXCJ0949.8+1708, the morphological analysis is not 
conclusive of the cluster dynamical status. We obtain w = 
0.0113 + 0.003, c = 0.197 + 0.007, and P3/P0 = (3.8+ 3.6) x 
10“^. As can be seen from Fig. [TTl re, and c are at the edge 
of the quadrant among radio quiet and radio loud clusters, 
and so does also P3/P0, within the error. We note that if the 
bias was taken into account in the P3/P0 computation, we 
would only be able to put an upper limit P3/P0 < 3.5x 10~^. 

The cluster sits in the central regions of the morpho¬ 
logical diagrams, in between regular and merging clusters. 
The values of w, c, P 3 /P 0 , and the presence of a BGC at the 
cluster centre, suggest that no major merger is occurring. 
However, present data do to allow us to exclude a major 
merger along the line of sight. Additional optical or weak 
leasing data are required to assess the cluster dynamical 
status. 

We have estimated the spectral index of the halo in 
RXCJ0949.8+1708 using NVSS and VLSSr ( VLA Low- 
Frequ ency Sky Survey Redux Source Catalog, iLane et ^ 
l2014| j images. In Fig. [8] the radio emission from NVSS is 
shown, superimposed onto the GMRT LR image. We have 
estimated the spectral index in the region where both im¬ 
ages have signal above three times the noise. Since it is not 
possible to estimate the contribution of the sources (A to 
H in Fig. [6} in the NVSS image, only an upper limit to 
the halo flux can be derived. Since no emission is present in 
the VLSSr image above three times the noise we have put 
an upper limit to the halo flux. Assuming that the spec¬ 
trum is described by a single power-law between 74 MHz 
and 1.4 GHz, we can constrain the spectral index in the 
range 0.6 < a < 1.5. Hence, we can only exclude a very 
steep spectrum radio halo. 


4.1.4 RXCJ 1354 . 6 -h 7715 

The results from the substructure analysis are P3/P0 = (4 + 
1) X 10“®, w = 0.080 ± 0.004, and c = 0.201 ± 0.006, which 
place the cluster among the merging ones, although the high 
value of c indicates that the core has not been disrupted by 
the merger. Alternatively, the high value of c could indicate 
a second merger along the line of sight. 

These result s are i n agreement with the optical analysis 
bv iHoresh et al.l (|2010h . 

No radio halo is detected at the sensitivity threshold of our 
observations, and under some assumptions (see Sec. [3] ) we 
cannot exclude the presence of diffuse radio emission a factor 
hve below the Pi. 4 GHz — SZ correlation. 


4.2 CIZAJ1938.3+5409: the first hadronic halo? 

The radio halo discovered in GIZAJ1938.3+5409 is the first 
direct detection of radio emission in a radio off-state cluster. 
As it is at least eight times under-luminous in radio for its 
SZ signal, we discuss two possibilities for the origin of the 
radio emission: 


(i) The observed emission in GIZAJ1938.3+5409 could 
be due to secondary electrons originating from hadronic 
i nteractions between thermal and non-therm al protons 
([Pfrommer X EnfilinI l2004bl : lEnfilin et al.l l201lh . We note 
that the mentioned works aimed at explaining radio 
halos with a higher power than the one observed in 
GIZAJ1938.3+5409. However, they can be easily extrapo¬ 
lated at lower radio powers because the particle spectrum 
predicted by the models is a simple power-law. The ex- 
tra polation at lower rad io powe r has been alredy done e.g. 
by iBrunetti X LazarianI (|201lh . According to this work, 
the electrons originating from hadronic collisions could 
produce radio emission that is a factor ~ ten below the 
Lx — Pi. 4 GHz correlation. Fermi upper limits indicate that 
the ratio of non-thermal G R protons to thermal e nergy 
densities is lower than few % (|Ackermann et giD 20j4L 


Sta cking a sample of clusters in radio. 


Brown et al.l 


(l201lh have claimed the first detection of the radio-off state 
of clusters. In order to better understand the possibility 
that this halo is produced by secondary electrons, we have 
investigated a set of simple models for GIZAJ 1938.3+5409. 
As a first approximation, we have assumed a /3-model for 
the gas density, rescaled to the mass of the cluster. We 
have assumed an isothermal temperature of T 6.4 keV as 
indicated by X-ray data, and we have computed the profile 
of cosmic ray (GR) protons by requiring that the energy in 
GR is a fixed fraction of the gas energy, ecu- 

Furthermore, we have modelled the magnetic f ield o f 
the cluster as we have done in iBonafede et al.l (|2013l i. 
Specifically, we have generated a power-law spectrum 
distribution of vector potential in the Fourier space for a 
200^ grid, randomly drawn from the Rayleigh distribution, 
and we have computed the magnetic field in the real space 
as B = V X A., ensuring V • B = 0 by construction. We 
have assumed that the maximum coherence scale of the 
magnetic field is 34 kpc, that the power-law of fluctuations 
is a Kolomogorv spectrum and that the average magnetic 
held stren gth in the central Mp c^ is (P) = 2 /xG, as found 
for Goma |Bonafede et ahllioiC 


Since the radio spectrum of the halo is not well con¬ 
strained, we explore a few possibilities for the spectrum of 
the cosmic ray protons, and adjusted the energy budget of 
GRs so that the secondary radio emission matches our obser- 
vation. We compute the ha d ronic y-ray em i ssion, following 
IPfrommer X EnfilinI (|2004al i ; IPonnert et al.l (|20101 ) , with the 
only difference that for the hadronic cross-section we use the 
parametrisat i on of the proto n-proton cross sect ion given by 
ikelner et al.l (|2006l L see also lVazza et al. (l2015ll for details . 
The secondary emission is based on Dolag X EnBlinI (|2000h . 

Results are listed in Table [3l All models with a spectrum 
of GRs p ^ —3 (with dN/dE oc E^) produce hadronic 
emission below the deepest limits by Fermi, which are 
at the level of ~ 10“^°ph/(cm^s) (e.g. Huber et al. 2013, 
Ackermann et al. 2014). The model with p = —3.1 produces 
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emission above the detectability of Fermi. However, the 
energy budget of CRs nearly equals the thermal energy, 
making the model unrealistic. 

A direct observation of this cluster with Fermi is un¬ 
practical due to its location in the Zone of Avoidance. Our 
calculations suggest that a hadronic origin for the radio 
halo in GIZAJ1938.3+5409 is consistent with the current 
Fermi limits. 

However, we note that if the cluster was at the luminosity 
distance of Coma (~ 99 Mpc) it should be detected in the 
7 -rays only if p ^ —2.3, i.e. a ^ 1.2. These estimates must 
be taken as approximative because of the uncertainties in 
both the magnetic held distribution and strength, which 
translate into an uncertainty in the normalisation of the 
energy budget of CRs for a given spectrum. For instance, 
for a magnetic held {B) =0.1 pG, Fermi upper-limits would 
be violated for p = — 2 . 1 . 

Our calculations suggest that if the halo has a hadronic 
origin, some clusters of the same mass as CIZAJ1938.3+5409 
in the nearby Universe might host radio and 7 -ray emission 
just above the current level of detectability. These halos do 
not represent the classic giant radio halo population, but 
rather a much broader population of “radio-quiet” objects, 
that still have to be detected. 

(ii) The second possibility to explain the radio emission 
in CIZAJ1938.3+5409 is that we are observing a radio halo 
in its rapid transition phase between the upper-limits re¬ 
gion and the Pi.4GHz-SZ correlation. A direct test, as done 
for the hadronic origin of the radio emission above, is not 
possible here because of the many free parameters. How¬ 
ever, we can make some statistical considerations in order 
to establish if the low power radio halo could be due to the 
same mechanism that powers radio halos on the Pi.4GHz-SZ 
correlation. 

In the framework of re-acceleration models, merger- 
induced turbulence lights up the radio emission, lifting the 
clusters from the radio off-state region of the Pi.4GHz-SZ 
plane up to the radio on-state, on the Pi.4GHz-SZ corre¬ 
lation. The steepest spectrum radio halos are expected to 
be below the correlatio n, as they represent t he latest stage 
of the radio halo life (jPonnert et al.l l2013l j . The halo in 
CIZAJ1938.3+5409 could be in the very early or very late 
stage of its radio-on phase. In both cases, its spectrum 
should be steep (a ^ 1.5). As no emission is detected in 
the VLSSr image, we can only put an upper limit to the 
spectral index a < 1.9. We cannot exclude that the radio 
halo has a very steep spectrum and that we are observing it 
just before it enters its radio-off state^ _ 

According to the simulations by IPonnert et al.l (|2013l j. 
clusters should spend ~ 0.1- 0.2 Gyr in the region of the 
Ri. 4 GHz — Lx plane between the upper limits and the radio 
power measured for GIZAJ1938.3+5409. Glust ers as massive 
as G IZAJ1938.3+5409 formed ~ 6 Gyr ago (|Giocoli et al.l 
bnOTfl . Hence, assuming that all the clusters follow the same 
path in the Pi .4 ghz — Lx plane, one would expect to find at 
maximum one cluster in 20-40 objects of the same mass 
and at the same redshift in the transition region, which is 
not at odds with current statistics. A proper estimate should 
account for the different merging rate at different redshifts. 
Hence, the number we hnd must be regarded as an upper 
limit. 


If the radio emission detected in GIZAJ1938.3+5409 is 
due to the same mechanism that powers the halo on the 
Pi. 4 GHz-SZ correlation, the cluster should be undergoing a 
merger, which has not left a strong imprint in the X-ray 
emission (see Sec. 14.11) . We note that the cluster lies in a 
rather empty quadrant of the morphological diagram c — w 
fFigfTTl). in agreement with the hypothesis that we are ob¬ 
serving it during a short phase, in transition between merg¬ 
ing and non-merging clusters. 

The two possible scenarios described above make very dif¬ 
ferent predictions: if we are observing a very rare object in 
its rapid transition between the radio quiet and radio-loud 
phase, we expect to observe few objects of this type in large 
cluster samples. Instead, if the halo has a hadronic origin, 
radio emission will be detected in every massive cluster (the 
upper-limits in Fig. HU). The detection of hadronic halos 
would also allow to constrain the energy budget in GR pro¬ 
ton in the IGM. 

Future deep radio sur veys performed e.g. with LOFAR 
(jvan Haarlem et ^l2013l j and later on with the Square K ilo- 
meter Array will shed light on this (jGassano et 


5 CONCLUSIONS 

In this work we have presented and analysed the radio emis¬ 
sion from four galaxy clusters selected from the Planck cat¬ 
alog (jPlanck Gollaboration et al]l201lh . The clusters are all 
massive, with M 500 > 6 x 10^^M© and are good candidates 
to host radio diffuse emission. We have used the same pro¬ 
cedure to reduce and analyse the radio emission from the 
clusters, and with the help of Chandra archival data we have 
performed a substructure analysis to quantify the dynami¬ 
cal status of the clusters. We have detected three radio halos 
and we have placed one new upper limit in the Pi .4 ghz — SZ 
plane. 

Detection of diffuse radio sources in three out of 
four newly observed clusters, or in four out of hve if we 
consider the previously published results on GL1821+643 
(Bo nafede et al. 2014) selected from the same SZ sam¬ 
ple (jPlanck Gollaboration et al.l l201lh . is indicative of a 
high fraction of occurrence of radio halos in SZ se¬ 
lected samples. Such high percentage have previously 
been suggested based on the analysis of NVSS radio 
images of Planck clusters (| Sommer k, Basul l2014l b The 
small sample size of the current work does not permit to 
make an y statistical predictions. However, the hrst Planck 
catalog (jPlanck Gollaboration et al] l201lh shows no net 
bias with respect to the later Planck SZ catalogs (e.g. 
IPlanck Gollaboration et ahl l2013l j] at high signal-to-noise 
(corresponding to M500 > 8 x 10^^M©. Hence, any con¬ 
clusion we may draw on the high radio halo fraction will be 
equally valid for later, more complete SZ catalogs. 

Our results can be summarised as follows: 

• We have discovered diffuse radio emission in three 
clusters, namely Abell 1443, GIZAJ1938.3+5409, and 
RXGJ0949.8+1708. In RXG J0949.8+1708 a tent ative radio 
halo was already claimed bv I Venturi et al.l (|2008l j 

• Glusters that fall in the same region of the morphologi¬ 
cal diagrams P3/P0, le and c can or not host radio halo, and 
its power can follow or not the correlation, indicating that 
the dynamical status is not the only cause of radio emission. 
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Table 3. Parameters of our models of hadronic secondary emission for PSZl G086.47+15.31, always assuming a Coma-like magnetic 
field configuration, with (B) = 2/rG in the central Mpc^. See text for details. 


model 

P 

a 

ecR 

-P 325 MHz 

[W/Hz] 

Po.5-100 GeV 
[erg/s\ 

Fo. 5-100 GeV 
[ph/{cm^s)] 

1 

-3.1 

1.7 

0.91 

2.4 X 1024 

5.2 X 1042 

3.9 X 10-10 

2 

-2.9 

1.6 

0.13 

2.4 X 1024 

8.5 X 1041 

6.5 X 10-11 

3 

-2.7 

1.5 

0.018 

2.4 X 1024 

1.4 X 1041 

1.1 X 10-11 

4 

-2.5 

1.3 

0.0093 

2.4 X 1024 

2.3 X 1040 

1.7 X 10-12 

5 

-2.3 

1.2 

0.0004 

2.4 X 1024 

3.6 X lO®*! 

2.7 X 10-10 

6 

-2.1 

1.1 

0.0001 

2.4 X 1024 

5.5 X 10“ 

4.2 X 10-14 


• The cluster Abell 1443 hosts peculiar emission which 
does not obviously fall into the categories of radio halo 
or radio relic. A bright P-shaped source is found at the 
cluster centre, which we interpret as a bright filamentary 
emission, simila r to t h e one found in MACSJ Q717+3745 
(jBonafede et al.l l20Q9l : Ivan Weeren et al.l l2009l b Alterna¬ 
tively, it could be a peculiar wide-angle tail radio galaxy. 

• The cluster CIZAJ1938.3+5409 hosts diffuse emission 
that we classify as radio halo. The total power of the halo at 
1.4 GHz, assuming a conservative a = 1.2 would be a factor 
eight below the Pi. 4 GHz—SZ correlations, placing the cluster 
in the “radio quiet” zone. The cluster could be caught in its 
rapid (~ 0.1 -0.2 Gyr) transition phase between the Pi. 4 GHz- 
SZ correlation and the upper-limit region of the Pi. 4 GHz-SZ 
plane, although no evident signs for an ongoing merger are 
visible in the X-ray image. 

Alternatively, it could be the first direct detection of a 
hadronic radio halo. We have verified that this hypothesis is 
consistent with the mos t recent 7 -ray upper limits by Fermi 
(jAckermann et al]l2014l b 

Information about the spectral index will be helpful to dis¬ 
criminate between the two possibilities or to point towards 
a different origin of the radio emission. It remains to be un¬ 
derstood, though, if this kind or low power radio emission is 
common to all massive clusters or not. In any case, this is 
the first direct detection of a cluster in its radio-off state. 

• The cluster RXGJ0949.8+1708 hosts diffuse emission 
on a Mpc scale that we classify as radio halo. The emission 
is elongated in the N-S direction and does not follow the X- 
ray emission from the gas. Although previous analysis based 
on optical and X-ray observations indicate that the cluster is 
not in a merging stage, the X-ray morphological estimators 
suggest that the cluster is borderline between merger and 
non merger clusters. 

• No diffuse radio emission is detected in 
RXGJ1354.6+7715 at the sensitivity reached by our 
observations. Gonsidering a halo size of ~ 700 kpc and a 
detection theshold of ‘^cflr we have placed an upper-limit 
to the radio emission which is a factor ~ 5 below the 
Pi.4GHz-SZ correlation. 
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